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I.  INTRODUCTION 


In  recent  years  there  has  been  considerable  emphasis  on  the  research  and 
development  of  various  electromagnetic-propulsion  devices.  In  addition  to 
technological  problems  associated  with  the  devices  themselves,  it  has  long 
been  recognized  that  the  electrical-power  requirements  necessary  for  driving 
these  accelerators  would  be  quite  severe.  In  particular,  it  is  necessary  to 
develop  a  power  source  which  is  reusable,  which  is  lightweight,  and  which  is 
capable  of  delivering  megajoules  of  energy  on  a  time  scale  of  a  few  millisec¬ 
onds.  Various  types  of  power  supplies  are  now  in  common  use,  but  all  have 
fairly  significant  disadvantages.  A  brief  summary  of  some  of  the  more  common 
sources  and  their  limitations  has  been  presented  elsewhere. 1 

One  type  of  power  supply  alluded  to  frequently  in  the  past,  but  which 
until  recently  has  not  received  much  detailed  study,  is  the  inverse  railgun 
(IRG).  The  basic  concept  upon  which  the  IRG  works  is  quite  straightforward. 

It  is  constructed  as  a  railgun  but  with  some  load  resistance  and/or  inductance 
connected  across  the  muzzle  and  with  a  seed  current  initially  in  the  circuit. 
Rather  than  the  current  providing  the  principal  force  that  drives  the  armature 
away  from  the  flux  entrained  in  the  bore,  however,  the  armature  is  driven 
toward  the  flux  by  an  external  force.  If  it  proceeds  rapidly  forward  the  cur¬ 
rent  rises  and  mechanical  energy  is  converted  into  electrical  energy.  Basic¬ 
ally  three  factors  must  be  present  for  the  IRG  to  operate  as  an  efficient 
electrical  generator:  First,  the  external  force  which  drives  the  armature 
must  be  sufficient  to  cause  current  gain  in  the  circuit.  If  the  velocity  is 
too  low  the  energy  will  be  dissipated  resistively  before  current  gain  can  oc¬ 
cur.  Second,  the  armature  must  be  completely  decelerated;  otherwise,  there 
will  be  inefficient  conversion  of  kinetic  to  electrical  energy.  Finally,  the 
electrical  energy  must  be  extracted  from  the  IRG  in  a  timely  manner  before  the 
armature  can  be  accelerated  backwards  by  the  magnetic  forces  in  the  circuit. 

In  fairly  recent  work,  Marshall^  has  considered  using  an  IRG  to  charge 
an  inductor  which  was  in  turn  to  be  used  for  earth-to-space  rail- launch  appli¬ 
cations.  He  concluded  that  the  device  was  worthy  of  further  study.  Later, 
Benford,  Kares,  Brooks  and  Goldman^  investigated  the  possibility  of  using  IRGs 
to  charge  inductive  stores  positioned  along  a  segmented  railgun.  At  least  for 
some  applications  IRGs  were  found  to  be  more  practical  for  this  use  than  the 
more  commonly  used  homopolar  generators.  Still  more  recently,  Powell  and  Jami- 
son^  considered  the  possibility  of  using  an  IRG  to  power,  first,  a  railgun 
and,  second,  a  load  of  constant  resistance  and  inductance.  (The  constant- 
impedance  load  might  be  used  to  model,  say,  a  plasma-heating  device  in  an 
electrothermal  gun. )  Whereas  the  latter  application  did  not  appear  to  be 
particularly  promising,  the  former  application  was  thought  to  merit  further 
study. 

In  this  report  we  will  be  concerned  with  extending  the  work  of  Ref.  3  to 
Investigate  in  detail  the  possibility  of  using  an  IRG  to  power  a  railgun  (RG). 
^Marshall,  R.A.,  "A  Reusable  Inverse  Railgun  Magnetic  Flux  Compression  Genera¬ 
tor  to  Suit  the  Earth-to-Space-Rail- Launcher,"  IEEE  Trans.  Magn.  MAG-20,  223 
(1964). 

^Benford,  J.,  Kares,  R.  Brooks,  A.,  and  Goldman,  E.,  "High  Energy  Railgun  De¬ 
signs,"  IEEE  Trans.  Magn.  MAG-20,  1*07  (1981*). 

^Powell,  J.D,  and  Jamison,  K.A.,  "Analysis  of  an  Inverse  Railgun  Power  Source," 
IEEE  Trans.  Magn.  MAG-22,  1669  (1986). 
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In  particular,  we  will  account  more  accurately  than  before  for  the  resistance 
of  both  sets  of  rails  and  both  armatures.  A  more  exact  treatment  was  consid¬ 
ered  necessary  since  our  previous  analysis  indicated  that  the  operation  of 
the  device  in  this  configuration  was  highly  sensitive  to  the  resistance  of 
both  the  IRG  and  the  load. 

The  report  is  organized  as  follows.  In  Sec.  II,  the  model  is  discussed 
and  the  governing  equations  presented.  In  Sec.  Ill,  an  exact  analytic  solu¬ 
tion  to  these  equations  is  obtained  for  the  case  in  which  resistive  effects 
are  neglected.  Most  of  the  material  in  this  section  is  presented  in  a  more 
abbreviated  version  in  Ref.  3.  In  Sec.  IV,  resistance  is  included  in  the  cal¬ 
culations.  Finally,  in  Sec.  V  the  results  are  summarized  and  some  indication 
of  future  work  is  given. 


II.  MODEL  AND  GOVERNING  EQUATIONS 

The  model  to  be  employed  in  the  calculations  is  shown  schematically  in 
the  upper  half  of  Fig.  1.  The  subscript  one  is  used  throughout  to  denote 
variables  associated  with  the  IRG  and  the  subscript  two  denotes  variables  as¬ 
sociated  with  the  RG.  A  seed  current  of  magnitude  i^  flows  initially  through 

the  two  armatures,  of  mass  m^  and  m^,  and  along  the  rails  in  the  direction 
indicated.  The  two  armatures  are  initially  located  at  x^  =  0  and  x^,  =  D, 
where  D  is  the  length  of  the  IRG.  At  time  t  =  0  an  external  force  is  applied 
to  m^  causing  it  tc  move  to  the  right.  If  the  motion  is  sufficiently  rapid 
the  current  in  the  circuit  will  rise  and  kinetic  energy  will  be  converted 
into  electrical  energy.  During  this  time,  m^  is  decelerated  by  the  J  X  $ 

force  (or  magnetic  pressure)  while  m2  is  accelerated.  The  object  is  to  deter¬ 
mine  for  what  parameters  the  velocity  v^  of  mg  can  be  maximized. 

A  more  conceptual  drawing  of  the  IRG-RG  configuration  is  shown  in  the 
lower  half  of  the  figure.  For  reasons  indicated  in  the  subsequent  analysis, 
the  most  desirable  design  is  usually  a  multi-turn  IRG  used  in  conjunction 
with  a  single-turn  RG.  Shown  in  the  drawing  is  a  three-turn  IRG.  Current  is 
conducted  along  the  rails  in  the  direction  indicated  by  the  arrows. 


It  is  convenient  in  all  subsequent  analysis  to  use  dimensionless  variables 
for  all  quantities  associated  with  the  operation  of  the  IRG.  The  manner  in 
which  these  quantities  is  made  dimensionless  is  shown  in  Table  I.  The  method 
of  normalization  was  rather  arbitrary  and  chosen  to  make  the  governing  equa¬ 
tions  as  simple  as  possible.  In  the  table,  L'  ls  the  inductance  per  unit 

length  of  the  IRG,  a  is  the  conductivity  of  the  rails,  u  is  the  free-space 
permeability,  and  v^  is  given  by 


v 


n 


2  E 


(2.1) 


where  Eg  is  the  electrical  energy  initially  in  the  circuit,  i.e.,  h  D* 
Physically,  v  is  the  velocity  reached  by  the  IRG  armature  if  it  was  started 
from  rest  and  accelerated  to  the  left  by  the  magnetic  force  of  the  seed  cur¬ 
rent,  with  no  energy  transferred  to  the  RG  armature. 
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TABLE  I.  Normalization  of  Variables 


Normalizing 

Factor 


Symbol  for  Dimen 
signless  Quant it 


Current 


IRG~ Armature  Position 


RG  Armature  Position 


IRG  Armature  Velocity 


RG  Armature  Velocity 


Resistance 


i:  v 
1  n 


Inductance 


IRG  Skin  Depth 


D/v 

n 

t  D 


HH' 

'■pav  J 
n 


RG  Skin  Depth 


wn 


Distance  Beneath  Rail 
Surface 


f — —I*5 


The  equations  which  govern  the  behavior  of  the  IRG-RG  configuration  fol¬ 
low  directly  from  the  equations  of  motion  of  each  armature  and  Kirchoff's  law 
for  the  circuit.  In  the  dimensionless  coordinates  just  defined  these  equations 
become 


=  vx  =  -I2/2  , 


(2.2) 


X2  "  V2  2  a  8  ’ 


(2.3) 


(L  I)  +  I  R  =  0 


(2.M 


The  dots  denote  differentiation  with  respect  to  the  dimensionless  time  x,  and 
L  and  R  denote  the  dimensionless  inductance  and  resistance  of  the  complete 
system  (see  Table  I).  It  is  evident  that  L  is  given  by 


.  j|«  *.v  klJlal*1Ai*  VJi*  tfltf 
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L  =  1  -  X1  *1  <*2  -  11 


(2.5) 


while  R  depends  upon  the  specific  form  of  the  resistance  assumed  for  the  IRG, 
the  RG,  and  the  two  armatures.  The  quantities  a  and  6  are  given  by 


=  m^/m^ 


(2.6) 


6  *  li/L2 


(2.7) 


where  L'  is  the  inductance  per  unit  length  of  RG.  In  the  most  general  case 
there  is  a  force  acting  upon  the  IRG  armature,  but  in  cases  to  be  studied  here 
it  is  assumed  for  simplicity  that  the  force  acts  impulsively  at  t  =  0,  provid¬ 
ing  m_  some  initial  velocity,  and  is  zero  afterwards.  The  force  is  therefore 
not  included  in  Eq.  (2.2). 

Equations  (2.2)-(2.M  must  be  solved  subject  to  the  initial  conditions 


(t  =  0)  =  0 


(2.8) 


(t  =  0)  =  /y 


x2  (t  =  0)  =  1 


v2  (t  =  0)  =  0 


I  (t  =  0)  =  1 


The  quantity  y  represents  the  ratio  of  the  initial  kinetic  energy  given  to  m^ 
and  the  initial  electrical  energy  associated  with  the  seed  current.  Thus,  we 
have 


y  =  v1q/L'  ig  D 


(2.9) 


where  v^  is  the  initial  (real)  velocity  of  the  IRG  armature.  As  will  be 
seen,  the  parameters  a, 8,  and  y  are  very  influential  in  determining  the  oper¬ 
ation  of  the  IRG-RG  system. 
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III.  SOLUTION  FOR  NO  RESISTANCE 


In  the  limiting  case  in  which  the  resistance  R  of  the  IRG-RG  configura¬ 
tion  is  neglected,  an  exact  analytic  solution  of  Eqs.  (2.2)-(2.M  can  be  ob¬ 
tained.  While  the  solution  is  not  useful  for  obtaining  quantitative  results, 
it  does  provide  some  insight  for  understanding  the  case  in  which  resistance 
is  included. 

To  solve  Eqs.  (2.2)-(2.h)  we  make  the  transformation 


X  =  x1  +  a  8 


Y  =  x2  -  6  x1  . 


(3.1) 


After  obtaining  a  first  integral  of  ( 2.h ),  we  find 


X  =  0  , 


Y=— 6(l  +  a62) 


(3.2) 


(3.3) 


Y  -  1  +  8 


(3.1*) 


The  solution  to  Eq.  (3.2)  obeying  the  initial  conditions  specified  in  (2.8) 
is  evidently 


X  =  /y  x  +  a  8 


Substituting  (3.1*)  into  Eq.  (3.3)  and  integrating  we  obtain 


(3.5) 


i  -  t  le’  y  *  (l - 6 — 

a  Y  -  1  +  8 


)]h  . 


(3.6) 


The  lower  sign  is  to  be  taken  in  Eq.  (3.6)  prior  to  the  time  Y  goes  to  zero, 
while  the  upper  sign  is  to  be  taken  for  later  times.  Finally,  integrating 
(3.6)  again  yields  r  versus  Y,  i.e.. 


t  =  ££  +  6(l/Vt0B2)  tanh_1(  8^7b)  ±  # 


.  e(l/a  +  82)  -1  pT 

+  rrr  tanh  /  r— 

b3/2  bv 


(3.T) 
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where 


b  =  62  Y  +  1/a  +  82  , 


u  =  82  y  (8  -  l)  -  1/a  -  8P  +  bY  , 


(3.8) 


(3.9) 


v  =  6  -  1  +  Y  . 


(3.10) 


One*3  X  and  Y  are  known,  x^  and  follow  from  hhe  transformation  inverse  to 
Eq.  (3.1),  namely. 


-  X  -  a  6  Y 

x,  =  - — 

1  1  +  a  62 


-  =  Y  B  X 


(3.11) 


1  +  a  B' 


(3.1?) 


It  is  evident  from  the  foregoing  analysis  that  the  operation  of  the  IRG- 
RG  system  depends  principally  upon  the  parameters  a,  8,  and  y.  One  is  normally 
interested  in  employing  the  device  for  values  of  a  which  are  fairly  small  and 
values  of  y  *  ich  are  fairly  large.  The  first  condition  is  necessary  to 
achieve  large  velocities  of  the  railgun  projectile,  and  the  second  condition 
is  necessary  to  achieve  significant  current  gain  in  the  system. 

It  is  furthermore  evident  that  for  a  given  a  and  y  the  parameter  8  can  be 
varied  so  as  to  optimize  the  velocity  imparted  to  the  projectile  m^.  It  is 
fairly  obvious  physically,  and  can  be  shown  mathematically,  that  the  optimum 
value  of  8  is  one  which  allows  the  IRG  armature  to  be  completely  decelerated 
just  at  the  time  it  reaches  the  end  of  the  IRG  gun  tube,  i.e.,  at  x^  =  D.  At 

later  times  its  velocity  is  negative  and  the  armature  accelerates  slowly  to 
the  left.  For  values  of  8  rfhich  are  smaller  than  the  optimum  value,  m^  will 
not  be  completely  decelerated  and  whatever  kinetic  energy  it  has  at  the  end 
of  the  gun  tube,  as  well  as  whatever  electrical  energy  is  stored  in  the  cir¬ 
cuit,  cannot  be  recovered.  On  the  other  hand,  for  values  of  8  larger  than 
optimum,  is  decelerated  too  fast  and  "bounces  back"  toward  the  left  with  a 

substantial  portion  of  its  kinetic  energy.  Again,  this  energy  cannot  be  re¬ 
covered.  In  practice,  8  can  be  varied  by  changing  the  configuration  of  the 
bore  and/or  the  number  of  turns  employed  by  one  or  both  guns.  In  subsequent 
analysis,  we  will  hold  L'  fixed  and  vary  L'  to  achieve  the  desired  value. 
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We  have  plotted  the  solution  represented  by  Eqs.(3.1l)  and  (3.12)  for  _ 
various  values  of  a  and  y  and  varied  8  until  the  maximum  (final)  value  of  v^, 

denoted  by  v^,  was  obtained.  Typical  of  the  results  are  those  shown  in  Figs. 
2-1*  in  which  are  plotted  the  dimensionless  variables,  x^,  x^»  v^,  v^,  and  I 
versus  the  dimensionless  time  t  for  three  different  sets  of  values  of  a  and 

Y* 

For  the  case  in  which  a  =  0.1  and  y  =  lU.O  (Fig.  2),  the  optimum  value  of 
8  was  found  to  be  about  3.8  with  both  higher  and  lower  values  of  6  producing 
smaller  final  velocities  of  m2*  As  can  be  seen  in  the  figure  the  final  veloc¬ 
ity  Vgf  is  about  11.7  or  about  three  times  the  initial  value  of  v^.  This 
velocity  is  reached  just  subsequent  to  the  time  m^  reverses  its  direction  at 
the  end  of  the  IRG  gun  tube_and  little  change  in  the  velocity  occurs  there¬ 
after.  The  final  value  of  X2  shown  in  the  figure  is  about  1*.7  so  it  is  evi¬ 
dently  necessary  that  the  railgun  be  nearly  four  times_  longer  than  the  IRG  to 
achieve  this  final  velocity.  Final  values  of  v^  and  x^  are  -0.69  and  0.86, 
respectively. 

Probably  the  most  significant  result  presented  in  Fig.  2  is  the  current 
profile.  The  small  width  of  the  pulse,  as  well  as  the  fact  that  the  maximum 
current  is  about  9*3  times  the  initial  value,  places  a  rather  severe  limita¬ 
tion  on  the  operation  of  the  IRG  for  these  values  of  the  parameters.  If,  for 
example,  we  choose  Ep  =  25  kJ,  m^  =  1  kg,  D  =  1  m,  and  L'  =  O.U  uH/m,  then 

v  =  2620  m/s,  i„  =  180  kA,  and  i  »  1.7  MA.  Such  large  currents  are 
2f  ’0  max 

likely  to  be  deleterious  to  the  rails  of  any  reasonably  sized  railgun.  Fur¬ 
thermore,  the  half-width  of  the  pulse  is  only  about  one-sixth  of  the  total 

acceleration  time  ('v  2.7  ms  for  the  real  parameters  Indicated  above).  Thus, 

the  force  on  the  armature  is  very  nonuniform  in  time.  This  situation  is  known 
from  conventional- gun  studies  to  be  very  undesirable;  not  only  is  damage  to 
the  projectile  likely  to  be  severe,  but  the  system  is  also  inefficient. 

The  high  maximum  current  can  be  considerably  reduced  provided  the  IRG 
armature  is  made  more  massive.  Such  an  effect  is  shown  in  Fig.  3  in  which 
are  plotted  results  for  the  same  value  of  y,  lk.0,  but  with  a  reduced  by  a 
factor  of  ten.  The  optimum  value  of  8  for  this  cas_e  turns  out  to  be  about 
12.0.  Final  velocities  and  positions  are  given  by  V2^.  =  36.9,  Xg^.  =  ^3.2, 

v^  =  -0.69,  and  x-^  =  0.85.  The  maximum  dimensionless  current  reached  is 

again  about  9.3  but  the  real  current  is  less  than  in  the  case  considered  pre¬ 
viously.  If,  for  example,  we  choose  Eg  =  25  kJ,  D  =  1  m,  and  =  O.b  uH/m, 

as  before,  but  choose  m^  =  10  kg  in  order  to  reduce  a  by  a  factor  of  ten,  we 

find  -  2608  m/'s,  i^  =  102  kA,  and  i^^  =  9^9  kA.  The  final  velocity  of 

the  projectile  is  therefore  about  the  same  as  bhe  previous  case  but  the  cur¬ 
rent  peak  has  been  reduced  by  nearly  a  factor  of  two.  The  half-width  of  the 
peak,  however,  is  still  small  relative  to  the  total  acceleration  time,  about 
8.7  ms  for  these  parameters. 


The  results  obtained  for  this  case,  compared  to  those  of  the  previous  ex¬ 
ample,  can  be  explained  approximately  as  follows.  Here,  the  more  slowly  mov¬ 
ing  projectile  will  compress  the  flux  entrained  in  the  bore  less  rapidly  than 
before.  Consequently,  comparable  current  gain  (1/1q)  and  a  retarding  force 
sufficient  to  decelerate  the  IRG  armature  can  be  achieved  only  if  L'  (or  B) 
is  increased.  Now,  however,  the  initial  current  1q  will  be  smaller  for  a 

fixed  stored  electrical  energy  and  the  maximum  real  current  will  be  corres¬ 
pondingly  smaller.  The  smaller  currents  also  lead  to  a  less  rapid  accelera¬ 
tion  of  the  projectile  m^  so  that  the  time  required  for  energy  transfer  to 

occur  is  much  longer  in  the  second  case.  Thus,  a  rather  long  RG  will  be  re¬ 
quired  and  the  additional  length  is  a  disadvantage  relative  to  the  previous 
case.  In  general,  a  compromise  between  very  high  currents  and  very  long  rails 
must  be  made.  This  point  is  discussed  further  when  resistive  effects  are  con¬ 
sidered. 

The  current  maximum  can  be  made  still  smaller  if  one  is  willing  to  employ 
smaller  values  of  the  energy  ratio,  y.  It  is  obvious  that  then  less  flux  com¬ 
pression  should  occur  and  the  current  gain,  i/ig,  should  be  smaller  than  for 

large  values  of  y.  Results  of  a  calculation  for  a  =  0.01,  as  before,  but  y 
reduced  by  50%,  i.e.,  y  =  7.0,  are  shown  in  Fig.  1*.  Here,  the  optimum  value 
of  B  wa_s  found  to  be_  about  10.8_and  final  values_of  the  pertinent  variables 
were:  =  27.0,  x^j.  =  17.5,  v^  =  -0.27,  and  xj^f  =  0.92.  The  maximum  rela¬ 

tive  current,  I  ,  is  now  only  about  It. 8,  i.e.,  about  a  factor  of  two  smaller 
„  max 

than  before. 

If  we  assume  =  10  kg,  D  =  1  m,  L'  =  O.k  pH/m,  as  before,  but  now  take 
Ee  =  50  kJ  so  that  y  is  half  the  value  considered  in  the  previous  example,  we 


find: 


v  =  2700  m/s,  i_  =  152  kA,  and  i  -  725  kA.  Thus,  the  maximum  cur- 
cl  u  max 


rent  has  been  reduced  by  about  25%  from  that  in  the  previous  case.  Again, 
however,  the  low  currents  require  exceedingly  long  RGs  to  effect  the  energy 
transfer.  Furthermore,  if  the  value  of  y  becomes  too  close  to  unity,  so  that 
the  initial  electrical  energy  is  comparable  to  the  initial  kinetic  energy  of 
the  armature,  the  IRG  becomes  a  rather  impractical  device;  the  entire  accel¬ 
eration  process  may  as  well  then  be  done  electrically. 

For  fixed  values  of  a  and  y  the  IRG  must  be  designed  so  that  the  optimum 
value  of  B  is  achieved.  It  is  therefore  important  to  know  how  sensitive  the 
final  velocity  v^f  is  to  the  variation  of  B.  To  partially  answer  this  ques¬ 
tion  we  have  plotted  in  Fig.  5  the  terminal  velocity  as  a  function  of  B. 
The  parameters  y  and  a  were  held  fixed  at  lit.O  and  0.1,  respectively,  as  for 
the  case  studied  in  Fig.  2.  For  values  of  6  to  the  left  of  the  vertical  line 
(6  -  3.8),  no  reversal  of  the  IRG  armature  occurs.  As  pointed  out  previously 
the  kinetic  energy  carried  by  m^  as  well  as  the  stored  electrical  energy  is 
then  lost  when  m^  reaches  the  end  of  the  gun  tube  and  the  system  is  very  in¬ 
efficient.  Consequently,  the  value  of  v ^  drops  quite  remarkably  as  one  pro¬ 
ceeds  to  the  left  of  the  vertical  line.  For  values  of  6  to  the  right  of  the 
line,  reversal  always  occurs.  As  B  increases  from  the  optimum  value,  however. 
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reversal  takes  place  farther  away  from  the  end  of  the  gun  tube  and  m^  acceler¬ 
ates  backwards  with  a  greater  portion  of  its  initial  kinetic  energy.  Again 
efficiency  is  lost.  Obviously,  however,  it  is  far  better  to  overestimate 
slightly,  rather  than  to  underestimate,  the  appropriate  value  of  B.  All  cal¬ 
culations  in  Fig.  5  were  terminated  when  =  1  if  no  reversal  occurred,  or 
when  v^j.  had  approached  a  terminal  value  when  reversal  did  occur. 

IV.  EFFECTS  OF  RESISTANCE 


A.  Analysis 

When  resistance  is  included  in  the  calculation,  the  situation  can  be  ex¬ 
pected  to  change  substantially.  In  particular,  for  a  given  value  of  L'  the 
current  gain  should  be  smaller  than  for  the  resistanceless  case  and  conse¬ 
quently  the  retarding  force  on  is  smaller.  It  will  therefore  be  necessary 
to  increase  the  inductance  ratio  B  in  order  to  decelerate  m^  to  zero  velocity 
just  at  the  end  of  the  IRG  gun  tube.  However,  the  same  basic  conclusions 
that  held  for  the  resistanceless  case  should  hold  here  as  well. 


We  consider  now  the  effects  of  resistance  on  the  performance  of  the  IRG 
configuration.  The  resistance  of  both  sets  of  rails  is  considered  as  well  as 
that  of  both  armatures.  It  is  assumed  that,  prior  to  time  t  =  0,  the  current 
in  the  IRG  is  brought  up  to  its  initial  value  of  iQ  by,  say,  discharging  a 

capacitor  in  a  time  tQ.  The  current  during  this  initial  discharge  time  is 
assumed  to  vary  sinusoidally  so  that  we  have 


i 


i0  sin 


TT  (t  +  tQ) 

77 


—  t  <  t 
x0  -x 


(fc.D 


and  m_^  and  m^  are  assumed  to  remain  fixed  at  x^  =  0  and  x^  =  D,  respectively. 
The  assumption  should  be  quite  reasonable  since  typically  t^  is  small  compared 
to  the  total  time  of  acceleration. 


Now  during  the  discharge  time  and  subsequently  after  the  armature  is  mov¬ 
ing,  the  resistance  of  the  IRG  rails  is  determined  by  the  IRG  skin  depth,  or 
the  extent  to  which  the  current  has  penetrated  beneath  the  rail  surface.  It 
is  furthermore  evident  that  while  m^  is  moving  forward,  the  skin  depth  is  in¬ 
dependent  of  position  along  the  rail  surface  and  the  total  resistance  is  given 

*y 


Ri  = 


2N2  (D  -  xx) 

"Vi 


(1..2) 


In  Eq.  (b.2),  N  represents  the  total  number  of  turns  in  the  IRG  (see  Fig.  l), 
<5^  is  the  skin  depth,  and  h^  is  the  total  height.  Thus,  each  turn  has  height 

h^/N  and  the  turns  are  series-connected.  Equation  (k.2)  does  not  hold,  of 
course,  once  the  armature  reverses  its  direction  and  the  skin  depth  now  depends 
upon  position  along  the  rails.  In  the  numerical  calculations,  however,  the 
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resistance  of  the  IRG  during  this  time  vill  be  neglected  since  it  is  generally 
small  relative  to  that  of  the  RG. 

It  is  shown  in  the  appendix  that  the  skin  depth  inside  any  of  the  conduc¬ 
tors  (IRG  rails,  RG  rails,  or  armature)  that  constitute  the  bore  of  the  IRG-RG 
system  is  given  at  any  time  by 


<5  =  - 


[n  fo]t 


(U.3) 


Here  n  denotes  a  unit  vector  normal  to  the  surface  (pointing  into  the  conduc¬ 
tor)*  B  is  the  magnetic  induction  field,  and  the  subscript  S  means  that  the 
quantity  in  question  is  to  be  evaluated  at  the  conductor  surface.  For,  say, 
the  upper  rail  in  the  IRG,  n  =  a  and  the  induction  field  for  y  >  0  is  ob¬ 
tained  by  solving  the  standard  diffusion  equation 


_3B  _  _1_  32B 
3t  uo  3y2 


(b.h) 


^tf fusion  in  directions  other  than  normal  to  the  conductor  surface  is  neg¬ 
lected.  The  boundary  condition  on  the  field  at  the  surface,  needed  to  solve 
the  equation,  is  taken  to  be 


Bs  =  B  (y  =  0)  =  uj2 


(^. 5) 


where  is  the  current  per  unit  height  on  the  rail  surface.  In  assuming  the 
above  results  we  have  neglected  any  dependence  of  the  induction  field  within 
the  bore  of  the  gun  on  z  and  taken  the  field  to  have  single  component  $  =  Baz. 

Such  an  assumption  is  not  expected  to  affect  greatly  the  skin-depth  calcula¬ 
tion.  If  we  now  write  Eqs.  (U.2)-(U.U)  in  terms  of  the  dimensionless  variables 
defined  in  Table  I,  we  find 


61  =  - m 


_3B  _  JFvB 

3t  "  3p 


CN2  h2  (1  -  x2) 

R  a  - — - 

1  hx  8  61 
In  Eq.  (1*.8) ,  C  is  given  by 


(b.6) 


(U.7) 


(U.8) 


L_  /  ub\H 

h2  V  Vjl  / 


(U.9) 


ft 

'A 

i*7il 

1 

s 

1 

$ 

si 

y, 

fe 

1 

£ 


B 

& 

1 

* 

|| 

h 


s 

$ 

<& 


| 

; 

& 


I 

i 

so. 


where  h2  is  the  railgun  height.  The  IRG  inductance  gradient  has  been  written 
as  8L '  since  we  customarily  vary  8  by  changing  L'  and  hold  as  well  as  the 
initial  electric  and  kinetic  energies  fixed. 

We  next  determine  how  the  inductance  gradient  L£  of  the  IRG  can  be  expect¬ 
ed  to  vary  as  one  varies  the  number  of  turns  and  changes  the  ratio 


si =  Vwi 


(I*. 10) 


where  Wj  represents  the  rail  separation  in  the  IRG.  As  Wq_  and  hj  are  varied 
it  is  assumed  that  the  cross-sectional  area  A^  of  the  IRG  remains  fixed.  For 
a  rectangular  cross-section,  single-turn  gun  the  inductance  gradient  can  be 
approximated  by  Batteh's  formula^ 


L1  =  G(sl} 
where  G(si)  is  given  by 

s  ^ 

G(Si)  =  s1  tan  1  (s1)  +  -jp  log  (l  +  1/s^)  -  log  (1  +  s2). 


(It. 11) 


(It. 12) 


For  an  N-turn  gun  of  the  same  geometry,  it  is  expected  that  the  inductance  L' 
would  become  multiplied  by  approximately  N8  in  keeping  with  standard  results 
for  square-bore  guns.  If  it  is  finally  assumed  that  the  RG  is  a  single-turn 
square-bore  gun,  then  we  obtain  from  Eq,  (It.ll)  and  the  arguments  above 


L'  ItN2  G(s  ) 

6  =  =  > 
L2  *  S1 


(It. 13) 


where  we  have  noted  that  G(l)  =  v/U.  Evidently,  then,  8  can  be  varied  by 
varying  the  turns  N  in  the  IRG  as  well  as  by  changing  the  height-to-separation 
ratio  s^,  until  the  optimum  value  is  found. 

The  skin  depth  of  the  RG  can  be  determined  via  considerations  analogous 
to  those  for  the  IRG.  Here  it  must  be  remembered,  however,  that  the  current 
has  not  been  flowing  the  same  length  of  time  at  all  points  along  the  RG  tube. 
In  particular,  the  current  is  initiated  at  any  point  only  when  the  projectile 
m^  arrives  at  the  point  in  question.  Consequently,  the  RG  skin  depth  6^  is  a 

function  of  x,  the  position  along  the  gun  tube,  and  the  total  RG  resistance 
must  be  found  by  integrating  each  contribution  along  the  gun-tube  surface. 

In  terms  of  the  usual  dimensionless  variables  we  have 

^Batteh,  Jad  H.,  "Momentum  Equation  for  Arc-Driven  Railguns,"  J.  Appl.  Phys. 
J56,  3182  (1981*).  See  also,  "Arc-Dynamic  Calculations  in  the  Railgun,"  Bal¬ 
listic  Research  Laboratory  Contract  Report  No.  ARBRL-CR-00521,  November 
1983. 
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The  skin  depth  6 ^  (x,t)  is  obtained  by  solving  (^.7)  at  various  points  x  along 
the  RG  gun  tube  and  using  Eq.  (U.6)  with  subscript  one  replaced  by  two. 


A  problem  appears  to  arise  in  Eq.  (l+.lk)  as  x  -*■  x^  since  <$2  0.  It  is 


possible  to  argue,  howeve£,  that  the  integral  in  C1*.^)  still  converges.  In 
fact,  in  the  vicinity  of  x^  the  skin  depth  is  well  approximated  by  the  expres¬ 
sion  valid  for  constant  current,  namely. 


$2  - 


(U.15) 


where  At  represents  the  time  the  current  has  been  flowing  at  the  point  in 
question.  If  this  expression  is  substituted  into  (it.ll*)  and  the  resulting  in¬ 
tegral  performed,  it  is  possible  to  show  that  the  contribution  to  the  rail 


resistance  between  some  point  x^  (near  x^)  and  X2  is  given  approximately  by 


tt.  2C  rV2  ^X2  ”  xm) 

R2  =  T  1 - * - J 


(1*.16) 


We  finally  consider  the  resistance  from  each  of  the  armatures  m^  and  it^. 

Here  it  is  assumed  for  simplicity  that  each  armature  has  the  same  conductivity 
as  the  rails.  It  is  apparent  then  that  the  current  diffuses  into  each  of  the 
armatures  at  the  same  rate  as  in  the  IRG  rails.  We  can  therefore  write 


2  Mi, 


—  (1  +  1/s, ) 


(1*.17) 


This  result  accounts  for  the  resistance  of  both  armatures,  but  neglects  resis¬ 
tive  effects  associated  with  the  rail-armature  Interfaces. 


The  equations  above  are  sufficient  to  yield  the  resistances  R^,  R2,  and 


R^  as  a  function  of  time,  and  the  total  resistance  R  in  Eq.  (2.1*)  is  simply 

the  sum  of  the  three.  Our  procedure  has_been  to_assume  a  value  for  B  and  solve 
numerically  the  resulting  equations  for  x^,  x2»  v^,  v 2,  and  I  using  the  assumed 

value.  The  parameter  B  was  then  varied  and  the  process  repeated  until  the 


maximum  value  of  Vg  was  reached. 
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Numerical  Technique 


In  the  numerical  solution  of  the  foregoing  equations,  a  simple  Euler 
method  was  employed  for  the  ordinary  differential  equations  (2.2)-(2.1),  and  a 
completely  explicit  technique  for  Eq.  (U.l).  The  time  step  dx  and  grid  spac¬ 
ing  dy  were  chosen  always  to  satisfy  the  stability  criterion5 


X  =  ■  4r—  <  1/2  .  (1+.18) 

(dy)2“ 

Typical  values  used  were  dx  =  0.0005  and  dy  =  0.05. 

A  dimensionless  length  of  the  entire  system,  denoted  by  ,  was  then 

max 

assumed  and  the  system  divided  into  200  segments  so  that  dx  in.  Eq.  (U.lU)  was 
given  by  3c  /200.  If  subsequent  calculations  revealed  that  x2  was  too 
max  max 

small  for  m^  to  reach  a  terminal  velocity,  or  if  it  was  significantly  too 

large,  the  value  was  changed  appropriately.  The  diffusion  equation,  Eq.  (^.7), 
was  then  solved  at  every  grid  point  a^Long  the  x  direction  in  the  RG  provided 
that  x  <  x  in  order  to  determine  62(x,  t).  We  normally  found  that  one  hund¬ 
red  grid  points  of  separation  dy  interior  to  the  rail  surface  were  sufficient 
to  prevent  the  field  from  diffusing  all  the  way  through  the  rails  during  the 
entire  calculation  time.  It  is  therefore  necessary  and  assumed  that  the  rail 
thickness  Is  sufficiently  great  to  prevent  complete joenetration  at  any  point. 

A  similar  calculation  was  carried  out  to  determine  6^  £or  the  IRG,  but  only  a 
single  point  along  the  x  direction  was  employed  since  61  is  independent  of  x.^. 

Once  the  skin  depths  6^  and  6^  have  been  calculated,  the  rail  resistances 
R  and  R2  follow  from  Eqs.  (b.8)  and  (k.lk).  The  integral  in  ( U.lh)  was  evalu¬ 
ated  by  using  the  trapezoidal  rule  for  all  grid  points  between  the  first  and 
the  one  just  prior  to  x2»  Between  the  last  grid  point  and  x2»  the  contribu¬ 
tion  to  the  resistance  was  approximated  via  (U,l6).  Finally  the  armature  re¬ 
sistance  was  found  by  use  of  Eq.  (1.17). 

The  above  procedure  of  solving  the  diffusion  equation  in  the  rails  and 
thereby  calculating  61  and  «2  was  carried  out  at  every  time  step.  The  result¬ 
ing  resistance  was  then  substituted  into  Eq.  (_2. 1)  and  the  solution  marched 
forward.  The  calculation  was  terminated  when  x^  became  equal  to  unity  (no 
direction  reversal),  or  when  v2  had  clearly  approached  some  asymptotic  value 
(after  reversal  had  occurred).  The  parameter  6  was  then  varied,  as  discussed 
before,  and  the  calculation  repeated  until  the  optimum  value  :of  B  was  found. 

^Carnahan,  B.,  Luther,  H.A.,  and  Wilkes,  J.D.,  Applied  Numerical  Methods 
(Wiley,  New  York,  1969),  Chap.  7. 
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C.  Results 

We  now  consider  the  results  of  some  calculations  undertaken  using  the 
above  procedure.  For  each  case  parameters  needed  to  be  specified  in  the  cal¬ 
culation,  as  well  as  some  of  the  more  significant  results,  are  indicated  in 
Table  II. 


TABLE  II.  Values  of  Parameters  in  Numerical  Calculations 

INPUT  PARAMETERS 


[see  Eq.  (i*.l)  ] 


Fig  6 

Fig.  8 

Fig.  9 

0.1 

0.01 

0.005 

ll*.0 

7.0 

7.0 

3 

1* 

5 

0.05 

0.05 

0.05 

2.0  cm 

2.0  cm 

2,0  cm 

50  cm2 

320  cm2 

320  cm2 

2.1*6 

1.81* 

2.19 

1  m 

0.25  m 

0.25  m 

RESULTS 

1.79 

1.09 

1.20 

9.5  cm 

19  cm 

20  cm 

6.3 

15.2 

22.5 

7.7 

18.3 

21*.  6 

2.8 

8.0 

10.0 

16.6 

5.1*0 

5.52 

Shown  in  Fig.  6  are  results  for  a  =  0.1  and  y  =  ll*.0.  (The  resistance¬ 
less  case  for  these  parameters  was  indicated  in  Fig.  2. )  The  cross-sectional 
area  A^  of  the  IRG  was  somewhat  arbitrarily  chosen  to  be  50  cm2  and  a  three- 

turn  IRG  was  found  to  produce  the  optimum  value  of  6  for  the  most  nearly 
square  bore  (s^  =  1.79).  The  optimum  value  of  the  inductance  ratio  was  found 
to  be  about  6.3,  nearly  a  factor  of  two  larger  than  in  the  resistanceless 
case.  As  indicated  previously,  the  higher  IRG  inductance  is  necessary  to 
decelerate  m^  when  resistive  losses  cause  the  current  to  drop. 


rvYi/Yvwj 
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It  is  of  some  interest  to  consider  the  results  of  Fig.  6  for  the  same 
reel  parameters  as  indicated  for  Fig.  2:  Eg  =  25  kJ,  =  1  kg,  D  =  1  m,  and 

L'  =  0.1*  uH/m.  We  then  obtain  v0„  ^  1730  m/s,  i„  =  1^0  kA,  and  i  ^  2.3  MA. 
There  is  therefore  a  reduction  in  projectile  velocity  by  about  35#  and  an  in¬ 
crease  in  the  maximum  current  by  about  35#.  The  decrease  in  velocity  can  be 
ascribed  to  the  loss  in  efficiency  because  of  resistive  dissipation,  and  the 
higher  current  to  the  larger  inductance  gradient  L'.  The  length  of  RG  needed, 
however,  is  only  about  1.8  m  (see  the  entry  in  Table  II)  as  opposed  to 

3.7  m  in  the  perfect-conductor  case.  It  is  rather  unlikely  that  2.0  cm-high 
rails  could  withstand  these  extremely  high  currents. 

It  is  next  of  interest  to  examine  how  the  potential  per  unit  length 
along  the  IRG  rails  varies  as  a  function  of  time.  We  have  therefore  plotted 
in  Fig.  7  the  quantity 


V  =  — - =- 

1  2  ( 1  -  X.. 


(1*.19) 


versus  t  for  the  calculation  shown  in  Fig.  6.  At  early  times,  when  the  total 
current  is  nearly  constant,  the  skin  depth  increases  roughly  as  /F  and  the 
potential  per  unit  length  decreases  slightly.  Alternatively,  one  could  say 
that  the  fairly  large  current  induced  at  the  rail  surface  when  the  induction 
field  is  first  turned  on,  decreases  as  the  field  penetrates  the  conductor  and 
reduces  the  field  gradient  at  the  boundary.  As  flux  compression  occurs,  how¬ 
ever,  and  the  total  current  rises,  J  at  the  surface  must  also  increase  as  a 
consequence  of  Lenz's  law  and  in  an  effort  to  exclude  the  field  'f'rom  the  con¬ 
ductor.  The  potential  therefore  rises  sharply.  Similarly,  once  the  maximum 

value  of  I  is  reached  and  4^-  becomes  negative,  a  current  density  in  the  reverse 

direction  is  induced  at  the  surface,  again  as  a  consequence  of  Lenz's  law. 

This  current  density  is  imposed  on  what  was  originally  there  but,  if  4^  is 

di 

sufficiently  large  in  magnitude,  the  current  density  at  the  surface  can  be 
oppositely  directed  from  that  which  initially  existed.  It  is  for  this  reason 

that  the  potential  becomes  negative  in  Fig.  7  for  negative  ~.  Evidently, 

then,  the  skin  depth  61  is  also  negative  as  is  the  IRG  resistance  R.^  The  un¬ 
usual  negative  resistance  results  from  the  uncommon  use  of  a  circuit  equation, 
with  a  constant  inductance  gradient  [see  Eq.  (2.b)],  to  describe  a  situation 
in  which  current  penetrates  beneath  the  surface  of  a  conductor.  In  reality, 
the  negative  resistance  accounts  for  the  decrease  in  flux  beneath  the  rail 
surface  as  the  current  drops. 

It  is  evident  from  the  results  of  Fig.  6  and  the  discussion  in  Sec.  Ill 
that  the  magnitude  of  the  current  pulse  needs  to  be  reduced,  and  that  such  an 
effect  can  be  accomplished  by  increasing  the  IRG  armature  mass  as  well  as 

decreasing  the  energy  ratio  y.  It  is  furthermore  obvious  that  decreasing  the 
resistance  will  lead  to  greater  efficiency  of  the  IRG.  The  resistance  can  be 
decreased  primarily  by  increasing  the  bore  size,  A^,  and  decreasing  the 

length,  D.  Results  of  a  calculation  in  which  these  changes  have  been  effected 
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are  shown  in  Fig.  8  with  pertinent  input  parameters  as  well  as  some  of  the 
more  significant  results  in  Table  II.  Specifically,  D  was  reduced  to  0.25  m 
and  increased  to  320  cm2.  Values  of  a  and  y  were  the  same  as  for  the  re¬ 
sistanceless  calculation  in  Fig.  It,  i.e.,  0.01  and  7.0,  respectively.  The 
four-turn  IRG  was  found  to  give  the  closest  approximation  to  a  square  bore 
(s^  =  1.09)  and  an  inductance-gradient  ratio  of  15.2  was  necessary  to  deceler¬ 
ate  m^. 

If  we  choose  as  real  parameters  m^  =  10  kg,  D  =  0.25  m,  Ee  =  50  kJ,  and 
L'  =  O.b  pH/m  (the  same  as  for  Fig.  h  except  that  D  has  been  reduced  as  dis¬ 
cussed),  then  v^  jv,  1830  m/s,  i^  257  kA,  and  lmnv  MA.  The  reduction 

in  velocity  (about  30%)  from  that  obtained  in  Fig.  4  results  from  a  decrease 
in  efficiency  because  of  resistive  losses.  The  higher  currents  result  from 
the  greater  flux  compression  associated  with  the  larger  inductance  ratio  B, 
as  well  as  the  higher  initial  current  caused  by  the  smaller  value  of  D.  The 
length  of  the  RG  needed  to  achieve  this  final  velocity  is  about  seven  times 
the  length  of  the  IRG  (see  Table  II ).  For  D  =  0.25  m,  this  length  is  1.75  m, 
a  reasonable  value. 

If  a  is  made  still  smaller,  say,  0.005,  smaller  peak  currents  can  again 
be  achieved  at  the  expense  of  longer  rails.  In  Fig.  9  are  results  of  a  cal¬ 
culation  for  y  =  7.0,  as  before,  but  a  =  0.005.  Additional  input  parameters 
are  the  same  as  in  the  previous  example,  except  that  a  five-turn  IRG  was  found 
to  produce  the  best  approximation  to  a  square  bore  (s^  =  1.20).  The  optimum 
value  of  6  was  found  to  be  6  =  22.5.  If  we  take  for  real  parameters  D  = 

0.25  m,  Ee  =  50  kJ,  and  L'  =  0.U  pH/m,  as  in  the  previous  example,  but  increase 
m^  to  20  kg,  then  we  find  v2f  a,  rfitO  m/s,  i^  2L  21°  hA,  and  imx  2L  1*2  MA. 

The  RG  length  needed  to  achieve  this  velocity  is  about  9D,  i.e.,  about  2.25  m 
for  D  =  0.25  m.  The  reduction  in  velocity  by  about  5#  from  that  obtained  pre¬ 
viously  can  be  attributed  to  the  higher  resistance  of  the  five-turn  gun. 


IV.  SUMMARY  AND  CONCLUSIONS 

It  appears  that  satisfactory  coupling  of  energy  between  the  IRG  and  RG  is 
theoretically  possible.  Other  types  of  load, 3  in  a  similar  switchless  con¬ 
figuration,  that  we  have  investigated  have  not  proved  nearly  so  satisfactory. 
The  principal  reason  that  the  IRG-RG  combination  works  reasonably  well  is  ap¬ 
parently  the  time  variation  of  the  RG  impedance.  In  particular,  both  the  re¬ 
sistance  and  inductance  are  very  low  at  early  times  when  flux  compression  is 
necessary  to  achieve  current  gain. 

The  major  disadvantage  of  using  the  IRG  to  power  a  railgun  appears  to  be 
the  high  peak-to-average  current.  This  problem  can  be  ameliorated  to  some 
extent  by  judicious  choice  of  the  initial  electrical  energy  in  the  IRG  and  a 
careful  choice  of  the  IRG  armature  mass.  Specifically,  small  values  of  both 
a  (ratio  of  RG  projectile  mass  and  IRG  armature  mass)  and  y  (ratio  of  initial 
kinetic  and  electrical  energies  associated  with  the  IRG)  produce  smaller-height 
larger-width  current  pulses. 


Other  possibilities,  not  investigated,  also  exist  for  flattening  the  cur¬ 
rent  profile:  For  example,  Marshall^  has  considered  a  configuration  in  which 
the  IRG  inductance  gradient  is  not  constant,  but  rather  is  larger  at  the 
breech  than  muzzle.  A  cursory  examination  indicates  that  this  design  would 
cause  the  current  to  rise  faster  and  hence  be  more  nearly  constant.  The  ef¬ 
fect  can  be  achieved  by  varying  the  number  of  turns  on  the  IRG  as  one  proceeds 
from  breech  to  muzzle.  It  is  also  possible  that  use  of  nonimpulsive  driving 
forces  could  produce  more  nearly  constant  current,  but  we  have  not  thus  far 
been  able  to  find  any  reasonable  force  profile  for  obtaining  the  desired  ef¬ 
fect. 

In  addition  to  problems  associated  with  the  current  profile,  other  areas 
appropriate  for  future  investigation  include  accounting  for  additional  losses 
and  considering  designs  which  allow  for  switching.  Among  the  losses  not  ac¬ 
counted  for  here,  for  example,  is  the  stray  inductance  between  the  IRG  and 
RG.  This  loss  was  not  included  in  our  analysis  since  no  reliable  estimate 
of  its  value  was  known.  A  number  of  possibilities  for  enhancing  energy 
coupling  between  the  IRG  and  RG  (and  other  types  of  load)  exist  if  switching 
is  permitted.  It  would,  for  instance,  be  of  interest  to  examine  the  effect 
of  switching  the  load  into  the  circuit  only  after  maximum  flux  compression  has 
been  achieved.  We  have  not  previously  considered  designs  which  allow  for 
switching  because  the  switching  itself  constitutes  a  major  practical  problem. 
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LIST  OF  MOST  COMMONLY  USED  SYMBOLS 


Symbol 


A1 

B 

D 

hl 

h2 

i 

*0 

I 

L/L 

L2 

"l 

m2 

N 

R/R 

VRi 

V*2 

Vra 


Vl/Vl 


V2/V2 


NOTATION:  IRG  =  Inverse  Railgun;  RG  =  Railgun 

_ Quantity _ 

IRG  cross  sectional  area 
Magnetic  induction  field 
Length  of  IRG 

Rail  height  (total)  of  IRG 
Rail  height  of  RG 
Current 

Initial  current  in  IRG 
Dimens ionless *  current 
Real/dimensionless  inductance 
IRG  inductance  gradient 
RG  inductance  gradient 
IRG  armature  mass 
RG  armature  mass 
Number  of  turns  for  IRG 
Real/dimensionless  total  resistance 
Real/dimensionless  IRG  rail  resistance 
Real/dimensionless  RG  rail  resistance 
Real/dimensionless  resistance  of  IRG  and  RG  armatures 
Ratio  of  IRG  rail  height  and  rail  separation 
Ratio  of  RG  rail  height  and  rail  separation 
Time 

Real/dimensionless  IRG  armature  velocity 
Real/dimensionless  RG  armature  velocity 


For  normalization  procedure,  see  Table  I 


a 


VI 


*1 


! 


a 


1 


V  /V 

V  1 


xi/xi 


x2/x2 


V5i 


VS2 


Normalization  velocity  [see  Eq.  (2.l)] 


Real/dimensionless  IRG  potential  per  unit  length  of  rail 


IRG  rail  separation 


RG  rail  separation 


Real/dimensionless  IRG  armature  position 


Real/dimensionless  RG  armature  position 


Real/dimensionless  distance  "beneath  rail  surface 


Ratio  of  RG  and  IRG  armature  masses 


Ratio  of  IRG  and  RG  inductance  gradients 


Real/dimensionless  IRG  skin  depth 


Real/dimensionless  RG  skin  depth 


Ratio  of  initial  kinetic  and  initial  electrical  energies 
in  IRG 


Conductivity 


Dimensionless  time 


V.v'.V.V 
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APPENDIX 

The  purpose  of  this  appendix  is  to  derive  Eq.  (2.4)  and  to  infer  from  the 
result  a  general  formula  for  the  skin  depth  6. 

Consider  Fig.  1  and  assume  as  discussed  previously  that  the  magnetic  in¬ 
duction  field  in  the  bore  of  either  gun  is  given  by 


B1 ,2  "  •*  h,2 


(A.l) 


where  j  represents  the  total  current  per  unit  height.  From  Maxwell's  equations 
one  also  has 


U  Y  P  3B 

V  X  E  =  -  — 


where  E  is  the  electric-field  intensity. 
(A. 2)  we  find 


/  e  • «  -  -  rsff 


(A. 2) 


If  Stokes'  theorem  is  applied  to 


(A. 3) 


The  integral  on  the  left  is  the  line  integral  over  some  closed  contour,  and 
that  on  the  right  is  over  the  surface  area  enclosed  by  the  contour.  Finally, 
from  Ohm's  law 


J  =  a  (E  +  v  x  B) 


(A. 4) 


where  v  represents  the  velocity  of  the  conductor  carrying^current  density  J. 
Thus,  along  the  rails  v  =  0,  but  along  the  two  armatures  v  is  given  by  the 
armature  velocity. 

We  now  substitute  Eq.  (A.U)  into  (A. 3)  and  apply  the  result  to  a  contour 
that  passes  along  the  inner  surface  of  each  armature  and  along  the  inner  sur¬ 
face  of  the  rails.  Evaluating  the  integral  in  the  counterclockwise  direction, 
we  easily  find 


j--dt-,1»1B1*,2  w2  b2  = 


-W^M1  <D-V 


-  U  ^  W2  (D  -  Xj)  . 


Equation  (A. 5)  can  be  written 


(A. 5) 
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I 

tV 

NS*. 
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I 
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where  L  is,  according  to  the  standard  definition,  given  by 

L  =  j  Js  B  •  dA  =  jH-  Wj  (D  -  Xj)  +  ^  W2  (x2  -  D)  .  (A. 7) 

It  is  then  apparent  that  Eq.  (A. 6)  is  just  Eq.  (2.U),  written  in  terms  of 
dimensional  variables,  provided 


(A. 8) 


We  now  use  the  additional  Maxwell  equation 


pJ  =  V  X  B  , 


(A. 9) 


and  the  elementary  expression  for  the  resistance  of  a  conductor  of  length 


1  dx 
R:o  l  h«(x) 


(A. 10) 


Here  h  represents  the  conductor  height  and  6  the  thickness  or  extent  to  which 
the  current  has  penetrated  beneath  the  surface.  It  is  assumed  that  6  may  vary 
with  position.  Equations  (A. 9)  and  (A. 10)  are  now  substituted  into  Eq.  (A. 8) 
and  the  integrals  performed  along  the  IRG  rails,  the  RG  rails,  and  each  of  the 
armatures.  If  the  resulting  expression  is  inspected,  it  becomes  readily  ap¬ 
parent  that  the  two  sides  of  (A. 8)  are  equal  if  6  is  given  on  each  surface  by 


<5 


[fa  •  n]s 


(A. 11) 


In  Eq.  (A. 11),  n  represents  the  unit  normal  to  the  surface  (pointing  into  the 
conductor),  and  the  subscript  S  denotes  that  the  quantity  is  to  be  evaluated 
at  the  surface  of  the  conductor  in  question. 


As  an  application  of  Eq.  (A. 11)  consider  the  simple  case  of  constant  cur¬ 
rent  and  calculate  6  for  the  upper  rail  in  the  IRG.  For  this  case  an  exact 
analytic  expression  exists  for  the  induction  field  inside  the  conductor  (y>0),6 
viz.. 


B(y,t) 


(A. 12) 


where  Erfc  is  the  complementary  error  function  given  by 

Powell,  J.D.,  "Thermal -Energy  Transfer  from  Arc  to  Rails  in  an  Arc-Driven 
Railgun,"  IEEE  Trans.  Magn.  MAG-20,  395  (1984).  See  also  Ballistic  Research 
Laboratory  Report  No.  ARBRL-TR-02530,  October  1983. 
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Erfc  (rn)  =  —  f  e'r  dr  . 

Sn  0 


(A. 13) 


If  we  now  note  that  n  =  a  ,  differentiate  (A. 12),  and  substitute  into  (A. 11) 
we  obtain  ^ 
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(A. 14) 


Equation  (A. 14)  is  the  standard  expression  for  constant-current  skin  depth. 
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